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Andreev tunneling through a double quantum-dot system coupled to a ferromagnet 
and a superconductor: effects of mean field electronic correlations 
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We study the transport properties of a hybrid nanostructure composed of a ferromagnet, two 
quantum dots, and a superconductor connected in series. By using the non-equilibrium Green's 
function approach, we have calculated the electric current, the differential conductance and the 
transmittance for energies within the superconductor gap. In this regime, the mechanism of charge 
transmission is the Andreev reflection, which allows for a control of the current through the ferro- 
magnet polarization. We have also included interdot and intradot interactions, and have analyzed 
their influence through a mean field approximation. In the presence of interactions, Coulomb block- 
ade tend to localized the electrons at the double-dot system, leading to an asymmetric pattern for 
the density of states at the dots, and thus reducing the transmission probability through the device. 
In particular, for non-zero polarization, the intradot interaction splits the spin degeneracy, reducing 
the maximum value of the current due to different spin-up and spin-down densities of states. Nega- 
tive differential conductance (NDC) appears for some regions of the voltage bias, as a result of the 
interplay of the Andreev scattering with electronic correlations. By applying a gate voltage at the 
dots, one can tune the effect, changing the voltage region where this novel phenomenon appears. 
This mechanism to control the current may be of importance in technological applications. 

PACS numbers: 73.23Hk, 73.63Kv, 74.45.+C, 74.78Na 



I. INTRODUCTION 

The interest in transport properties of mesoscopic sys- 
tems has increased a lot due to their potential for present 
and future technologies. Recent advances in the experi- 
mental development of nanostructures are mainly aimed 
at the study of purely quantum phenomena and effects 
based on electron-spin properties (spintronics) . In par- 
ticular, hybrid resonant structures composed by one or 
more quantum dots (QD) coupled to normal (N), ferro- 
magnetic (F) and superconductor (S) metals have been 
studiedi*2iMiSi2iZiS. In systems composed by one quan- 
tum dot, electron-spin properties have been extensively 
explored. In the special case of junctions composed by a 
ferromagnet and a superconductor it is possible to con- 
struct spin valves which control the current flow through 
those systems. Andreev reflection permits such control, 
by varying the polarization of the ferromagnet attached 
to the system, as shown in several papers ^ 10 ' 11 ! 12 ' 13 ! 14 ' 15 . 
Andreev reflection^ (AR) is a mechanism in which a 
Cooper pair is formed in the superconductor from the 
combination of an incident electron coming from the nor- 
mal metal with energy uj and spin <r, with another elec- 
tron with energy — uj and spin —a. Both electrons enter 
the superconductor as a Cooper pair, leaving a reflecting 
hole in the ferromagnetic electrode. Andreev states are 
located within the superconductor gap, where no quasi- 
particles states are available. 

In this work we have studied the transport properties 
of a hybrid nanostructure composed by a ferromagnet, 
two quantum dot a 14 i 17 ' 18 i 19 i 20 i 21 i 22 , and a superconductor 
connected in series (F—QD a — QDb — S). The addition of 
an extra quantum dot will allow us to study the interplay 
of electron correlations at the dots (for both, intra and in- 



terdot interactions), with the Andreev current. Figure [T] 
shows a schematic diagram of the system. The supercon- 
ductor chemical potential is fixed to zero (us = 0) and 
the bias is applied to the ferromagnetic electrode. There 
are also applied gate voltages at the dots a and b, namely 
V ga and V g b, respectively. By using the non-equilibrium 
Green's function 9 - ! 23 ! 24 ' 25 , we have calculated the current 
(I), differential conductance (dl/dV), Andreev transmit- 
tance (Tar) and the local density of states (LDOS) at the 
dots. All quantities are calculated for energies within the 
superconducting gap, the relevant range for the Andreev 
reflection, as functions of the voltage bias. We have also 
included intradot and interdot Coulomb correlations at 
the dots, and have analyzed its influence on the electric 
current through a mean field approximation. In solids, 
both correlations compete to form charge or spin modu- 
lated structures. Those symmetry broken states are not 
possible in finite systems, as it is the case in our double 
dot sample 2 ^. However, dot a is coupled to a ferromag- 
net, which breaks spin symmetry, and dot b is coupled to 
a superconductor, which acts as a charge reservoir. Thus, 
interesting effects are expected, when the electronic in- 
teractions at the dots are taken into account. In this 
paper, those effects are displayed by the differential con- 
ductance, which shows asymmetric regions of negative 
values as a function of the applied bias, when the I x V 
characteristics are obtained. Negative differential con- 
ductance (NDC) have been observed in hybrid nanos- 
tructures composed by normal metals^ 7 -, semiconductor 
based devices 2 ^ and more recently in molecular Joseph- 
son junctions 29 . There are also some theoretical studies 
on the NDC effect in those systems, using models be- 
yond the mean field approximation^iMLi^iHiMi^ii^. In 
our work, electron interactions at the dots are treated 
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within a mean field approach. This approximation, plus 
additional correlations introduced through couplings to 
the F/S electrodes, gives rise to NDC effects. For An- 
dreev currents, correlation parameters at the dots are 
limited by the size of the superconductor gap. 




FIG. 1: (Color Online) Schematic diagram showing the sys- 
tem studied in this work. The dot coupled to the ferromagnet 
electrode (F) is called a, and b is the one coupled to the super- 
conductor (S). The superconductor has its chemical potential 
fixed to zero, and the voltage bias is applied to the ferromag- 
net. Gate voltages are also applied at the dots. The different 
couplings are also indicated in the figure. 

This paper is organized as follows: In section If 
we present the model under consideration and derive 
the transport properties by using the non-equilibrium 
Green's functions. In section III the numerical results 
are presented and discussed. Some conclusions are given 
in section IV. 



II. MODEL AND FORMULAE 

A. Hamiltonian 

The system displayed in figure Q] is described by the 
following Hamiltonian: 



H = Hf + 7~(s + 7~td q d + Ht, 



(1) 



TLp is the Hamiltonian of the ferromagnet F described 
by the Stoner model. The spin bands of F are shifted by 
h, the exchange energy. The ferromagnet chemical po- 
tential is fixed by the applied bias, i.e., [ip — eV. Tis is 
the Hamiltonian for a BCS superconductor with chemi- 
cal potential fixed to zero as the ground, [is = 0. Ttdqd is 
the Hamiltonian for the quantum dots in the mean field 
approximation, which permits an exact equation for the 
self-energy. The energies E aa and Eba are renormalized 
by the interactions JC (interdot) and U (intradot). The 
interactions also couple the renormalized energy levels 
with the mean occupations (h a ) and (fib)- In addition, it 
is included a gate voltage at the quantum dots a and 6, 
namely V ga and V g b , respectively. Ht is the Hamiltonian 
which describes all the tunneling processes: between dot 
a and the ferromagnet, with amplitude tf, between dots 
with amplitude t a b, and between dot b and the supercon- 
ductor, with amplitude t s . 



B. Green's functions 

To calculate the transport properties we have used the 
non-equilibrium Green's function methoc ffi 23 ' 24 i 25 . All 
the physical quantities can be cast in terms of the Green's 
function of the dots. By using the Nambu 4x4 notation 
the retarded Green's functions of the quantum dots are 
given by: 
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H-dqd = E a(J n aa + 2^ Fbafiba 



where, 



E aa = eVga + — {n b ) + — {n aa ) 

_ T . K, ,„ , U , 
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In these equations G^ a is the Green's function of the 
quantum dot a; G£ b is the Green's function of the quan- 
tum dot b; g aa and g£ b are the Green's functions of 
the dots a and b isolated from the electrodes; t a & de- 
scribes the coupling between the dots; SJ, and Eg are 
the retarded self-energies which describe the coupling of 
the dots with the superconductor and ferromagnet elec- 
trodes, respectively. Explicitly these self-energies are 
written as, 
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with r fa = 2ir\tf\ 2 Na is the coupling strength, with tf 
being the tunneling amplitude and N a the density of 
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states for the ferromagnet spin a band; and 
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where T s = 27r|£ s | 2 iV s (0), with iV s being the density 

of states of the superconductor in the normal state 

and p s is the modified BCS density of states p(oj) = 

U0(M-A) uO(A-\u\) . , , . 

— , with the imaginary part 



A 2 



^VA 2 



accounting for Andreev states within the gap2i2!. 

Besides the retarded and advanced Green's functions, 
it is necessary to obtain the Keldysh Green's functions, 
which are calculated by the equation of motion technique. 
Since it is used a mean field approximation for the inter- 
action, the result for this Green's function is exact. The 
equation obtained for the Keldysh Green's function of 
dot a is given by: 



with the "lesser" self-energy 



(10) 



£<H = £<M +ti 6 G£E|( W )G£(a;)t o6 . (11) 

Correspondingly, the Keldysh equation for quantum 
dot b is given by: 



where fs = f(co) is the Fermi function for the supercon- 

is the conventional 



ductor electrode and p = 

Vw 2 - A 2 
BCS density of states. 

Equation (TIT]) shows that the dot a, which is coupled 
to the ferromagnetic electrode on its left side, 'sees' on its 
right side an effective electrode as a result of the interplay 
of dot b with the superconductor. Equation (fT3|) can be 
interpreted in similar terms for dot 6, with a 'bare' su- 
perconductor electrode on the right side, and an effective 
electrode on the left, resulting from the interaction of dot 
a with the ferromagnet. Since the superconductor and 
the ferromagnet present different band structures, there 
is an intrinsic asymmetry in this system which manifests 
itself in the transport properties. 



C. Physical Quantities 

The Green's functions of the last section, calculated by 
the equation of motion method, permit to determine all 
the physical quantities necessary to analyze the transport 
properties of the F — QD a — QDb — S system. Since the 
interaction couples the dot levels through the mean occu- 
pation, as shown by equations ([2]) and Q , it is necessary 
to perform a self-consistent calculation to determine the 
occupation at the dots first. Then, one can proceed to 
calculate the physical quantities of interest. 

In the following we show the expressions we have used 
to compute the LDOS, the current, the transmittance 
and the mean occupation. 



G< (u) = G£ b HS< h HG£ b (a, 
with the "lesser" self-energy £^ b : 



(12) 
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The "lesser" self-energy for the ferromagnet electrode 
is given by: 
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in which f F = f(ui — eV) and f F = (ui + eV) are the 
Fermi functions for electrons and holes, respectively. 

The "lesser" self-energy for the superconductor elec- 
trode is given by: 



(15) 
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1. Local density of states (LDOS) 

The LDOS of the quantum dots comes from the ma- 
trix elements [11] and [33] of the retarded Green's func- 
tion matrix (electron components in Nambu space). The 
LDOS for dots a and b are, respectively: 



LDOS-A = — Im[GL,n + G^. 33 ] (16) 
LDOS-B = --lm[Gl b<11 + Gfc i3a ]. (17) 

7T 



2. Transmittance and current 

Since the current is conserved, it can be calculated at 
any point of the circuit. Here, we choose to calculate the 
current at the ferromagnetic electrode, as the temporal 
variation of the number of electrons , i.e.: 



dt 



where N F — J2ka a \a a k<j- By using the Heisenberg equa- 
tion and the definition of the "lesser" Green's function of 
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the dot a, it's possible to write the current as follow: 



III. RESULTS AND DISCUSSION 



I=^Jdcu [GL»£<H + G<»££H + H.C.] u+33 

(18) 

where the index 11 + 33 indicates a sum over the electron 
components in the Nambu space matrix. By substitut- 
ing the matrix elements, the current can be cast to the 
following form: 



- / dw A(u)(f F 



Ff). 



(19) 



In this work we only consider Andreev transport, 
for energies within the superconductor gap. Thus, the 
current amplitude corresponding to the contribution of 
quasi-particles tunneling is zero. 

The expression for the amplitude A(oj) is given by: 



Next, we present the results obtained from numerical 
calculations. Firstly, we show the local density of states 
(LDOS) at the quantum dots in the absence of electronic 
correlations. We investigate the effects of the different 
couplings of the model, namely the coupling between dots 
and the coupling of the dots with the electrodes. In the 
following, we discuss the role of the interdot interaction. 
In the specific case with no polarization in the ferromag- 
net, P = 0, we have observed the appearance of the NDC 
effect, for some values of the interaction. The inclusion of 
the intradot interaction lifts the spin degeneracy, as can 
be seen by the splitting of the differential conductance 
peaks. In this case the NDC also appears, but is reduced 
with the increase of the polarization. 



A. Noninteracting case 



^ = r /T (|GL, 1 4l 2 r/T + |GLa2l 2 r 



fi) 



r /i (l G aa,34| 2r /T + \ G aa, 32^ 



fl) 



The transmittance is obtained from the current for- 
mula: 

Tar = i[r /T (|G; Q , 14 | 2 r /T + |GL, 12 | 2 r n ) + (20) 

r /l (l G oa,34| 2r /T + l G aa,32| 2r /l)]- 



3. Self- consistent calculations 

Since the Green's functions are dependent on the mean 
occupations via equations © and (j3J), it is necessary to 
calculate those quantities at the dots. From the definition 
of the "lesser" Green's function, one straightforwardly 
obtains the system of equations below: 
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These integral equations have to be solved numerically 
in a self-consistent way. Once the occupation numbers 
are obtained, it is possible to calculate the other physical 
quantities. Results are shown below. 



In this section, the LDOS of the quantum dots is de- 
scribed without interactions. This permits us to ana- 
lyze the resonance structure presented by these quanti- 
ties, which plays a central role in all transport proper- 
ties, including the interacting case as well. The system 
is asymmetric, since the ferromagnetic left electrode is 
modeled with a continuous density of states, while the 
superconductor electrode, on the right side, presents a 
gap for quasi-particles states, with a complex density of 
states within the gap, corresponding to evanescent An- 
dreev states, responsible for the Copper pair conversion 
at the interface. 

Figure shows the LDOS at the dots, for different 
values of the interdot coupling (electron hopping be- 
tween dots) tab, in units of the superconducting gap. 
For t a b = 0.02, the dots are almost decoupled from each 
other. As a result, features of the density of states mainly 
reflect the coupling with the electrodes. LDOS for dot a 
presents one peak centered in to = with a finite width. 
The broadening results from the hybridization of the dot 
level with the ferromagnetic band. There is a finite prob- 
ability for the electron to escape from the dot to the elec- 
trode. On the other hand, LDOS-B presents two sharp 
symmetrical peaks. This resonant structure represents 
the hybridization between the dot level with the Andreev 
states. The peaks corresponds to the electron and hole 
channels, as expected from the BCS model for the super- 
conductor electrode. 

As the coupling between the dots is increased, two ad- 
ditional peaks at the center emerge for both LDOS, as 
observed in the examples for t a b — 0.30 and 0.46. These 
peaks come from the resonance between the discrete dot 
levels. For LDOS-B, the intensity of the Andreev peaks 
decay with the dot coupling. 

The effect of the coupling with the superconductor 
is illustrated by the LDOS curves shown in figure [3] 
When the interaction with the superconductor is weak 
r s = 0.05, both LDOS present a two-peaks structure 
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(Color Online) LDOS for different values of coupling between dots, t a t- Fixed parameters: Ff — 0.1, T s — 1.0, and 

0.02, the LDOS 



k B T = 0.01 and eV ga 



Cl- 



ub 



0. (a) LDOS for QD a, coupled to the ferromagnet. For t a b 



is dominated by the coupling with the ferromagnet, revealed by the broadening of the central peak. By increasing t a b, the 
resonances from the superconductor (external peaks) and from interdot coupling (central peaks) appear, (b) LDOS for QD b, 
coupled to the superconductor. For t a b = 0.02 the LDOS is dominated by the coupling with the superconductor, as displayed 
by the equidistant peaks around u) = 0; these peaks are the Andreev resonances. By increasing t a t, the interdot coupling peaks 
(central peaks) appear in addition to the Andreev resonances. All parameters are expressed in units of the superconductor gap. 
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FIG. 3: (Color Online) LDOS for different values of the coupling with the superconductor, F s . Fixed parameters: T/ = 0.1, 
tat = 0.5, P — 0, ksT = 0.01 and eV ga = eV go = 0. (a) LDOS for quantum dot a, coupled to the ferromagnet. (b) LDOS for 
quantum dot 6, coupled to the superconductor. The distinct behavior of LDOS-A and LDOS-B is explained by the stronger 
coupling of b with S. The increase in the coupling with S results in a bigger separation between the resonance peaks. All 
parameters are expressed in units of the superconductor gap. 
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FIG. 4: (Color Online) LDOS for different values of coupling with the ferromagnet, Tf. Fixed parameters: T s = 1.0, tab = 0.5, 
P = 0, feT = 0.01 and eV ga — eV g b = 0. (a) LDOS for quantum dot a, coupled to the ferromagnet. (b) LDOS for quantum 
dot b, coupled to superconductor. By increasing Tf, the density of states presents a broader pattern, displaying an admixture 
between the ferromagnetic energy band and the hybridized states of the dots. All parameters are expressed in units of the 
superconductor gap. 



resulting from the interdot coupling (t a f, = 0.5 in the 
examples shown). When the coupling with the super- 
conductor is increased, the Andreev peaks appear and 
are more intense in LDOS-B. 

The effect on the LDOS by varying the coupling with 
ferromagnetic electrode is shown in figure |4j When the 
coupling with the ferromagnet is increased, the discrete 
structure of the LDOS is transformed into a continuum 
of states, as a result of the hybridization of the discrete 
dot levels with the continuous band of the ferromagnet. 
Internal peaks almost disappear for Tf > 0.5. 

Isolated quantum dots present one level degenerate in 
spin. When coupled to each other with t a b, there is an 
admixture of them, resulting in a bonding and an anti- 
bonding levels, in analogy with a H2 molecule^. In our 
model, those levels corresponds to the central peaks of 
the LDOS. When the electrodes are attached to the dou- 
ble dot system, the above peaks broaden and two addi- 
tional peaks appear corresponding to the superconduct- 
ing Andreev states. By tuning the parameters of the 
model, it is possible to change the number of peaks, their 
widths and the distance between them, which in turn can 
be used to control the current. 



B. Interacting case: Inter-dot interaction 

In figure[5ja) we plot some Ix V characteristics, for dif- 
ferent values of the interdot interaction /C. These curves 
show a plateau pattern which is due to the peak structure 
of the LDOS at the quantum dots. When the interaction 



is increased, the plateau value is reduced, ranging from 
/ = 0.90 for K = to I = 0.30 for K = 0.45, since higher 
values of the interaction implies a stronger Coulomb re- 
pulsion between dots. But for small voltages [eV < 0.30), 
we observe an unusual behavior, where the trend is in- 
verted, although this is a tiny effect. In figure EJb), we 
plot the corresponding differential conductance, which 
allow for a better resolution of the I x V curves. The 
symmetric structure for K, = is broken when K, 7^ 0, 
the asymmetry being more pronounced the higher the 
values of JC. For some examples of the figure, NDC in 
the characteristics is found around JC — 0.6. From our 
numerical calculations, NDC effects are present in the 
range 0.08 < K < 0.4. For K greater than 0.4, NDC is 
suppressed and a positive peak emerges in dl/dV, as can 
be seen in the example for JC = 0.45. 

From these results, we conclude that the mechanism of 
the NDC is not linearly related to the Coulomb blockade 
effect. The interaction plays a more subtle role in chang- 
ing the transmittance of the system. In fact, looking at 
the differential conductance, we note that when increas- 
ing }C, the second peak for positive bias is suppressed. 
Thus, for some values of fC there is a suppression of some 
of the resonant peaks, and this causes an additional re- 
duction of the transmittance for some values of the ap- 
plied bias. This effect causes the differential conductance 
to assume negative values. The Andreev transmittance 
is displayed in figures [5tc) and (d) . There is a variation 
with the applied bias, in contrast to the non-interacting 
case. The interaction couples the occupation number at 
the dots, which implies a non trivial dependance of the 
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FIG. 5: (Color Online) (a) current (7) versus applied bias (V) for some values of interdot interaction, (b) Corresponding 
differential conductance, showing regions of negative values, (c) Andreev transmittance (Tar) for some values of applied 
bias for tC = 0.22. (d) Andreev transmittance (Tar) for some values of /C, for applied bias eV = 0.99. Fixed parameters: 
Tf = 0.19, F s = 0.40, t ab — 0.5, P = 0, U — 0, ksT = 0.01 and eV ga = eV g t = 0. All parameters are expressed in units of the 
superconductor gap. 
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FIG. 6: (Color Online) LDOS for different values of interdot interaction. Fixed parameters: F/ = 0.19, F s = 0.40, t a t = 0.5, 
P = 0, U = 0, ksT — 0.01, eV = +0.62 and eV ga = eV g b = 0. (a) LDOS for quantum dot a, coupled to ferromagnet. (b) 
LDOS for quantum dot b, coupled to superconductor. In both, the interaction introduces an asymmetry which is related to 
the NDC effect. All parameters are expressed in units of the superconductor gap. 
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transmittance with the applied bias. In figure [5^c), we 
plot the transmittance for JC = 0.22, for some values of 
the applied bias. There is a reduction of the amplitude 
with the increase of the bias, but the spectrum is sym- 
metric with respect to w. In figure [5](d) , we show the 
transmittance at fixed bias eV = 0.99, for various values 
of JC. There is a reduction of the transmittance and a 
shift of the peaks, however the variation is not system- 
atic, as shown by the example for JC = 0.45, which does 
not follow the trend of the other values. Reduction of the 
transmittance with increasing bias is one of the causes of 
NDC. However, the absence of NDC for negative bias 
strongly hints that there are additional ingredients to 
explain the phenomenon. One important factor is the 
asymmetry in the LDOS, which appears when the inter- 
action JC is turned on. In figures a) and (b), we show 
the effect of the interaction on the peak structure of the 
LDOS at the dots. The symmetric shape seen in figures 
[U [21 and H] is lost when the interaction is included. By 
increasing the interaction towards the gap value, some 
peaks are suppressed (one central and one external) and 
other are reinforced (one central and one external), the 
LDOS presenting a more localized character. Central 
peaks of the LDOS are associated to states resonating 
between dots, while external peaks are the channels for 
the Andreev reflection. The symmetry is critical to al- 
low electronic transfer through the structure, since the 
sum of the energies of the electrons available to form a 
Cooper pair have to be equal to the chemical potential of 
the superconductor (which is zero). Thus, the Andreev 
current is optimized when the LDOS peaks are symmet- 
ric, and the suppression of one of them causes an effective 
reduction of the current, with the emergence of the NDC 
effect. 

When a negative bias is applied to the ferromagnetic 
electrode, its chemical potential is reduced in relation to 
the superconductor's one. Thus, the current is estab- 
lished by extracting Cooper pairs from superconductor 
electrode. Those electrons, with antiparallel spins, fill 
the lower energy states available at the dots, so only the 
peaks of the negative frequency branch of the LDOS will 
participate in the conduction process. This is the expla- 
nation for the absence of NDC for negative bias in the 
I x V characteristics, as shown in figure a). 

Next, to show that the NDC effect originates from 
the asymmetry of the LDOS, we have recalculated the 
JC = 0.22 case of figure [5l but now applying a gate volt- 
age at dot a, while keeping the gate voltage at the other 
dot fixed and equal to /is = 0. The results are qualita- 
tively similar if the gate potential at dot b is varied, while 
the one at a is kept fixed and equal to /xs- As shown in 
figuresEIa) and (b), the NDC appears for negative values 
of the applied bias, for V ga approximately ranging from 
-0.26 to -0.39. In the range -0.13 < V ga < +0.13, the 
NDC appears for positive bias. To make contact with 
the asymmetry of the LDOS at the dots, in figures [5](a) 
and (b), we plot the LDOS-B for values of the bias at the 
threshold of the NDC, namely eV = -0.61 and +0.62. 
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FIG. 7: (Color Online) Current and differential conductance 
vs. applied bias for some values of the gate voltage at dot a. 
Fixed parameters: Tf = 0.2, T s = 0.4, t ab = 0.50, P = 0, 
JC = 0.22, U = 0, k B T = 0.01 and V gb = 0. (a) current versus 
applied bias: the gate potential modifies the current profile 
appearing some regions of NDC in both, negative and positive 
values of the applied bias, (b) Differential conductance for the 
corresponding I x V curves. All parameters are expressed in 
units of the superconductor gap. 



The LDOS-A presents a similar behavior. By tuning the 
gate voltage V ga , we can change the amplitude and posi- 
tion of the peaks in the LDOS. Eventually, some of the 
peaks vanish, given rise to the NDC effect. The above 
figures corroborate the role of the LDOS in the appear- 
ance of the NDC regions. In fact, peaks of the LDOS and 
transmittance are resonances resulting from the coupling 
between the dots and the electrodes. As shown in figures 
[5]and[ni the interdot interaction affects the LDOS and the 
transmittance in a way similar to "destructive interfer- 
ence" , changing the position and amplitude of the peaks. 
One can tune such "destructive interference" , by intro- 
ducing a gate voltage, thus modifying the values of bias 
where NDC takes place, as shown in figure [5) This pro- 
cess of controlling the current through the device, may 
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co (A) 

FIG. 8: (Color Online) LDOS of dot b vs. electron energy for 
some values of the gate voltage at dot a. Fixed parameters: 
V f = 0.2, T s = 0.4, t ab = 0.50, P = 0, k B T = 0.01, V gb = 0, 
K, = 0.22 and U = 0. (a) LDOS-B for applied bias eV = 
—0.61. The second peak for negative energy is progressively 
suppressed and vanishes at eV ga ~ —0.39. (b) LDOS-B for 
applied bias eV = +0.62. The first peak for positive energy, 
that was absent for eV ga = 0, emerges with application of 
the gate voltage. All parameters are expressed in units of the 
superconductor gap. 



be important for practical applications. 



C. Interacting case: Intradot interaction 

Finally, we present mean field results that include the 
intradot (onsite) interaction hi, with no intcrdot repul- 
sion (K. = 0). As shown in equations ([2]) and ([3]), the 
intradot interaction splits the up and down-spin states 
at each quantum dot, with the corresponding splitting of 
the transmittance and differential conductance peaks. 

However, as indicated by equations ^ and ©, this ef- 
fect can only be observed for different up and down-spin 



i 1 1 1 1 1 1 1 1 1 1 1 r 
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FIG. 9: (Color Online) Current and differential conductance 
vs. applied bias for some values of the intradot interaction. 
Fixed parameters: F f = 0.1, F s = 1.00, t ab = 0.50, P = 0.80, 
K, = 0, k B T = 0.01 and eV ga = eV gb = 0. (a) Current 
vs. applied bias, (b) Differential conductance. The intradot 
interaction lifts the spin degeneracy, producing a splitting of 
the peaks in the dl /dV curves. Very different values of V j and 
r„ have been used in the example, to get a larger separation 
between the resonance peaks. All parameters are expressed 
in units of the superconductor gap. 



occupations. This condition is met for non-zero values 
of the polarization P of the ferromagnet, when different 
numbers of spin up and spin down electrons are injected 
to the dots. The rates at which electrons are injected are 
T n /h = T f (l + P)/h and T n /h = T f (l~P)/h, for spin 
up and down respectively. In figure[5]the I xV character- 
istics and the corresponding differential conductance are 
shown for different values of the intradot interaction. As 
long as W ^ 0, the peaks start to split, and for U = 0.90 
the differential conductance presents a clear pattern with 
eight peaks. The I x V characteristic, for U — 0.90, also 
shows a number of additional steps and a final plateau 
with a reduced value of the current. The reduction of 
the maximum value of the current with P is explained 
by the reduction of the available conducting channels, as 
discussed in reference [38l Since the current is established 
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by Andreev reflection, it is necessary an equal number of 
spin-up and spin-down electrons to form Cooper pairs. 
Since the density of states for spin down is smaller, the 
current is limited by the number of spin-down electrons. 

In the examples presented in figure O which corre- 
sponds to P = 0.90, NDC effects are absent. By increas- 
ing the polarization the NDC region is reduced, and even- 
tually disappears, when we further increase the polariza- 
tion. The mechanism that accounts for the NDC effect 
for intradot interaction is the same as the one presented 
in the previous sections: the reduction of the transmit- 
tance with the applied bias, combined with asymmetries 
of the LDOS. For values of the polarization of f» 0.90, 
the mean number of electrons participating in the con- 
duction is so reduced, that a further reduction of the 
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FIG. 10: (Color Online) Current and differential conductance 
vs. applied bias for different values of the ferromagnet po- 

Vgb 



larization. Solid lines: 



V ga = -0.39 and V gb = 



Vga = 

-0.10. 



-0.13. Dotted lines: 



Fixed parameters: Vf — 0.20, 



r s = 0.26, t ab = 0.40, K = 0.25, U = 0.25 and k B T = 0.01. 
(a) Current vs. applied bias, (b) Differential conductance. 
The increase of the polarization suppresses NDC by reducing 
the available states in the conduction process. All parameters 
are expressed in units of the superconductor gap. 



channels does not imply in a reduction of the electrical 
current. This is the cause of the absence of NDC in the 
examples of the figure 



D. Interaction case: 
Intradot and Interdot interaction 



In the last sections we have presented NDC results with 
one of the interactions absent. However, it is possible to 
observe NDC when both interactions are active. In the 
figure [TU] is shown some I x V curves for U = JC = 0.25. 
In spite of the high polarization values, P = 0.80 and 
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FIG. 11: (Color Online) Corresponding LDOS for some values 
of the applied bias. By adjusting the hopping t a b an d the 
coupling with superconductor r s , it is possible to reduce the 
peaks of the LDOS allowing to observe the NDC even at high 
values of the polarization. Fixed parameters: P = 0.80, Vf = 
0.20, V s = 0.26, t ab = 0.40, K. = 0.25, U = 0.25, k B T = 0.01 
and Vga — Vgb = —0.13. (a) LDOS for dot a, coupled to 
ferromagnet. (b) LDOS for dot b, coupled to superconductor. 
All parameters are expressed in units of the superconductor 
gap. 
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0.90, there are regions of NDC in all these curves. For 
e Vga = eVgb = —0.13 (solid curves), NDC appears for ap- 
plied bias eV > 0.48. For V ga = -0.39 and V gb = -0.10 
the NDC appears for negative bias eV < —0.48 (dotted 
curves). Besides the interactions and gate voltage values, 
another difference from the figure [9] are the values of the 
coupling constants t a b and T s . In fact, by changing the 
coupling constants, it is possible to make some peaks of 
the LDOS so small that the current can be sensitive to 
the reduction of the channels even in polarization values 
close to unity and the NDC can be recovered. In fact, in 
the corresponding cases of the figure [Til the LDOS dis- 
plays some peaks almost totally suppressed. An example 
is shown in the figure [TT1 where is plotted the LDOS for 
P = 0.80 and eV a = eVb = —0.13. When we increase 
the applied bias there is a suppression of the first and 
third peaks localized at u> = —0.54 and lu — 0.23, re- 
spectively. The suppression is almost complete for the 
first peak of the LDOS-A and for the third peak in the 
LDOS-B curve for eV = 0.57 and eV = 0.83. Since the 
Andreev reflection requires a symmetric pair of channels 
in order to conduct, the process is dominated by those 
suppressed peaks allowing to observe the NDC in cases 
with high values of P. 

IV. CONCLUSION 

In this work, we have studied the effects of the in- 
tcrdot and intradot interactions on the transport prop- 
erties of double quantum-dot system coupled to a fer- 
romagnet and a superconductor. Energy parameters of 
the theory are limited by the size of the superconductor 
gap. This way, the conduction through the device is con- 
trolled by Andreev scattering processes. In the first part 
of the paper, the role of the coupling between dots and 
the coupling of dots with the electrodes was elucidated. 
Next, we study the effects of electronic correlations at 
the dots, within a mean field approximation. For both 
interactions, inter and intra-dot correlations, we found 
regions of negative differential conductance (NDC). Cor- 



relations tend to localized the electrons at the double-dot 
system, changing the LDOS at the dots by suppressing 
some peaks and shifting their positions, leading to an 
asymmetric pattern for the LDOS. 

This asymmetry reduces the number of available states 
to conduct through the Andreev reflection mechanism. 
The above phenomena, combined with the transmittance 
reduction with the applied bias, produce the NDC effect 
for some regions of the voltage bias. By applying a gate 
voltage, one can tune the effect and change the bias re- 
gion where NDC appears. Such kind of devices, as the 
one considered here, are on the verge of being produced 
by present technology, and our theoretical study may be 
useful to control the current in practical applications. K, 
and hi are intrinsic parameters, which are sample depen- 
dent. However, their effect can be monitored by the gate 
voltages, as shown in this contribution. With the addi- 
tion of a second ferromagnetic electrode, one may open 
the possibilities of crossed Andreev reflections and con- 
trol of the current by the relative polarization directions 
of the two ferromagnets. 

The results presented in this work were obtained from 
a mean field theory and have to be interpreted with at- 
tention. Indeed, one could argue that the NDC could be 
washed out by fluctuations. However, the results shown 
in the figures [10] and [IT] were obtained for high values of 
polarization and nonzero gate voltages where the mean 
field approximation works well since the high values of 
the polarization reduce the spin fluctuations. Therefore, 
we believe that the NDC is a real effect and not only a re- 
sult of the approximation. However, the exact extension 
of the validity of the approximation used in this work can 
be addressed only by experiments. 
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